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ABSTRACT 

We present a detailed study of the Thallium, neutron, Mercury and deuteron electric dipole 
moments (EDMs) in the CP-violating Minimal Supersymmetric extension of the Standard 
Model (MSSM). We take into account the complete set of one-loop graphs, the dominant 
Higgs-mediated two-loop diagrams, the complete CP-odd dimension-six Weinberg operator 
and the Higgs-mediated four-fermion operators. We improve upon earlier calculations by 
including the resummation effects due to CP-violating Higgs-boson mixing and to threshold 
corrections to the Yukawa couplings of all up- and down-type quarks and charged leptons. 
As an application of our study, we analyse the EDM constraints on the CPX, trimixing 
and Maximally CP- and Minimally Flavour- Violating (MCPMFV) scenarios. Cancellations 
may occur among the CP-violating contributions to the three measured EDMs arising 
from the 6 CP-violating phases in the MCPMFV scenario, leaving open the possibility of 
relatively large contributions to other CP-violating observables. The analytic expressions 
for the EDMs are implemented in an updated version of the code CPsuperH2 . 0. 
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1 Introduction 



With the imminent advent of the LHC, we are entering an exciting era for probing new 
physics at the TeV scale. If new physics is indeed observed at this scale, the questions of 
its flavour and CP structure will immediately become very critical. The non-observation of 
the Thallium ( 205 T1) [1], neutron (n) [2], and Mercury ( 199 Hg) [3] electric-dipole moments 
(EDMs) already provide remarkably tight bounds on possible new CP-violating phases 
beyond the Cabibbo-Kobayashi-Maskawa (CKM) one of the Standard Model (SM). Com- 
plementary to the direct explorations at the LHC, a new generation of precision low-energy 
experiments is also expected to play an important role. The new precision experiments will 
place much stronger indirect constraints on the possible CP and flavour structure of models 
of TeV-scale physics. In particular, if the proposed experiment searching for a deuteron 
( 2 H + ) EDM achieves the projected sensitivity [4,5], this will improve the existing bounds 
on possible CP- violating chromoelectric operators by several orders of magnitude [6] . 

One of the theoretically best-motivated scenarios of new physics is Supersymmetry 
(SUSY) [7]. Its minimal realization, the Minimal Supersymmetric extension of the Standard 
Model (MSSM), with SUSY broken softly at the TeV scale, addresses the naturalness of the 
gauge hierarchy, predicts gauge-coupling unification, provides a viable candidate for Cold 
Dark Matter (CDM) and may help explain the baryon asymmetry in the Universe (BAU) 
via a first-order electroweak phase transition [8] . An essential role in the generation of the 
BAU could be played by the new CP-odd phases that appear in the MSSM [8]. However, 
the non-observation of EDMs severely constrains these CP- violating phases [9] . 

The aim of this paper is to present a detailed study of the 205 T1, n, 199 Hg and 2 H + 
EDMs in the CP-violating MSSM. We include in our study the complete set of one-loop 
graphs [10,11], the dominant Higgs-mediated two-loop diagrams [12-14] of the Barr-Zee 
type [15] and the Higgs-mediated four-fermion operators [13,16,17], originally studied by 
Barr [18] within a two-Higgs doublet model E. We improve upon earlier calculations [19,20] 
by calculating the complete CP-odd dimension-six Weinberg operator [21] and by including 
resummation effects due to CP-violating Higgs-boson mixing [22] and to threshold correc- 
tions to the Yukawa couplings of the up- and down-type quarks and charged leptons [23]. 
We then use this compendium to derive representative constraints on the CP-violating 
parameters of phenomenologically relevant benchmarks in the MSSM, such as the CPX 
scenario [24], the trimixing scenario [14,25] and the general Maximally CP- and Min- 
imally Flavour- Violating (MCPMFV) framework [26], for selected sets of values of the 
CP-conserving parameters. Clearly, sufficiently small values of the CP-violating param- 
eters must be compatible with the experimental upper limits on the EDMs but, as we 
shall illustrate with explicit examples, larger values may also be allowed by non-trivial 

*However, we do not include here the contributions of a possible non-zero CP-violating QCD vacuum 
parameter 9. 
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cancellations. 

For the presentation of our analytic results, we follow the conventions and nota- 
tions of CPsuperH [27], especially for the masses and mixing matrices of the neutral Higgs 
bosons and SUSY particles. We note parenthetically that the new version of CPsuperH, 
CPsuperH2 . 0, includes an improved treatment of Higgs-boson propagators and Higgs cou- 
plings, and enables numerical predictions for a number of flavour-changing-neutral-current 
(FCNC) .B-meson observables, including CP-violating effects. On the basis of the results 
in this work, we further improve the code CPsuperH2 . by implementing the computation 
of the Thallium, neutron, Mercury and deuteron EDMs in the CP-violating MSSM. 

The layout of the paper is as follows. Section 2 presents all formulae relevant to the 
one-loop contributions to the electric and chromoelectric dipole moments of the charged 
leptons and quarks that result from chargino-, neutralino-, and gluino- mediated diagrams. 
Non-holomorphic threshold effects on the light-quark Yukawa couplings have been appropri- 
ately resummed, as these are the dominant source of higher-order corrections. In Section 3, 
we calculate the CP-odd dimension-six three-gluon Weinberg operator, taking into account 
loop diagrams involving t and b quarks and Higgs bosons, in addition to the previously 
known loop effects due to t and b squarks and gluinos. In addition, we present analytic 
results for the Higgs-mediated four-fermion operators and the dominant Higgs-mediated 
two-loop diagrams. In Section 4 we compute the 205 T1, n, 199 Hg and 2 H + EDMs in the 
CP-violating MSSM. Section 5 presents illustrative constraints on key soft SUSY-breaking 
parameters and CP-odd phases in the CPX, the trimixing and the MCPMFV scenarios. 
We summarize our conclusions in Section 6. 

2 One-Loop EDMs of Leptons and Quarks 

At the one-loop level, the charged leptons, e, \x and r, and the light quarks, u, d and s, can 
have EDMs induced by charginos, neutralinos and gluinos. The u, d and s quarks may also 
develop chromoelectric dipole moments (CEDMs) via the corresponding squark and gluino 
loop diagrams. In this Section we exhibit analytical formulae for the one-loop EDMs of 
charged leptons and light quarks, and the CEDMs of the quarks. 

We denote the EDM of a fermion by d^ and the CEDM of a quark by . The 
relevant (C)EDM interaction Lagrangian is given by 



where F^ v and G aiiv are the electromagnetic and strong field strengths, respectively, and 
the T a = X a /2 are the generators of the SU(3)c group. The interaction Lagrangian (12.11) 
leads to a matrix element of the form: 




l -dfF^fa, ul5 f - l -d c ;G a ^ 



q (y^ihT a q , 



(2.1) 



iM 



df e(q) ■ (p + p) u(p) 7s u(p) , 



(2.2) 
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Figure 1: Generic Feynman diagram for the EDM {dj) x of the fermion f induced by x 
exchange effects. The photon line can be attached to the sfermion f line, or to the x line 

ifX = xt- 



where p = q + p' and e(q) ■ p = e(q) ■ p', since e(q) ■ q = 0. We use the convention 
= | [7^, 7"] = i(Yl u -9^)- 

To set our coupling notations and normalisations, we write down the generic inter- 
action of a chargino Xi^ neutralino X? 994 or ghuno g a , collectively denoted by x-, with a 
fermion / and sfermion f \ 2 , as follows!!! 



C 



xff 



9 x l{( U. /'/./) /;• + gffi- (\ ;/',,/•) I'' + h.c. 



xff 



;2.3) 



Likewise, the interaction Lagrangians for the couplings of a photon A^ 1 with x an d /' read: 



C 



-eQ x (xinx)Ar 



Cf'f'A 



-ieQ f ,f* d, f'A» 



(2.4) 



Employing ( 12.31) and (12.4ft and taking into consideration the diagrams of Fig. [Tj we calculate 
the one-loop fermion EDM, 



dfy 



where 



rn 



167r 2 m|, 



xff'Wxff 



9Rij 91 



■ 13 



Q x A{m 2 Jm),) + Q f B(m 2 Xi /mj f ) 



(2.5) 



A(r) 



2 1 -r 



2 lnr 

1 — r 



B(r) 



2r lnr 



2(1 -r 



(2.6) 



with A(l) = —1/3 and B(l) = 1/6. We have checked that our analytic expressions for the 
one-loop EDMs are in agreement with [9] and [11]. 



^Here the convention for the couplings ql and ga is different from that used in [9]: Ql = R*k and 
9R = L* 
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We now present the individual one-loop contributions of charginos to the EDMs 
of charged leptons (df/e)*, up-type quarks (d^fe)* and down- type quarks (df/e)* . In 
detail, these are given by 



'd\ 



d^V* 1 ^™v± „ c± 



E^ 9m [(^; d )*^ W 1 [Q^A{mym\) + Q~ d B{my m \)\ ,(2.8) 



\ e J IQir 2 mi 

(f) x = ^e^^k^ 1 )*^ 

where the electric-charge assignments for the loop particles are: Q^± = ±1, = 2/3, 
= —1/3, and 

giT = -g(C R )iu glT = h*(C L ) i2 , (2.10) 

g£f = -g(c L Uu% + h d (c L )Uu%, g£f = K(c R )Uu% , (2.11) 
giT = -g(c R )a(u% + h u (c R uu%, = K(c L ) i2 (u% . (2.12) 

We note that the coupling coefficients defined in fl2.10p - fl2.12p appear in fl2.7p - fl2.9l) . re- 
spectively. 

Correspondingly, the contributions of neutralinos to the EDMs of charged leptons 
(df/e)*°, up-type quarks (d^/e)^° and down-type quarks (df/e)*° may conveniently be 
expressed as 

(f)" = ^E^^^V^lO/^/m^), (2.13) 

with f — l,u, d. The neutralino-fermion-sfermion couplings are 

g?// = -V2gT(N: 2 {U%-^/2gt w {Q f -T()N: i {U%-h f N^ 

gf!/ = V2gt w Q f N a (U%-h}N ia (U%, (2.14) 

where the Higgsino index a = 3 (/ = /, d) or 4 (/ = u), T^' = — 1/2 and T% = +1/2. 

In addition to charginos and neutralinos, gluinos also contribute to the quark EDMs 
(df'/e)^. In the gluino mass basis, (df/e)^ is given by 

I ' d E \ 9 1 I m I 

[fj =^E^^[(^)*^^(WK), (2-15) 



where the gluino-quark-squark couplings are given by 

^e-^(U%, gf] - 



gil = -^e-^ 2 ^, gf] = +^=e + ^ 2 (U%. (2.16) 



As well as the EDMs, the chargino, neutralino and gluino loops can produce non- 
vanishing CEDMs for the quarks, (d^) x± ' x °^ . Their individual contributions are as follows: 





9s \ - 

167T 2 ^ 

hi 




9s 

16tt 2 4^ 

h3 


a q=u,d ) 


9s \ - 

16tt 2 ft 


(dq= u ,d) 9 = 


9s \ - 

8vr 2 V 



dj 

*i C*„,IY„ X <?<A* „X 991 D/' 2 /_2 \ 



9j< 



where 



°f)-i(r^(*-» + T^-T^)- < 218 > 



so that (7(1) = 19/18. 

Finally, it is important to stress that there are non-holomorphic threshold corrections 
to the light-quark Yukawa couplings that occur in the chargino and neutralino couplings. 
These corrections are proportional to the strong coupling a s , and become significant at large 
tan j3. We therefore resum these effects by redefining the light-quark Yukawa couplings as 
follows: 

y/2m u 1 , \plm c 1 



VSp l + Au/tp' VSp l + Ac/fy' 

1 + A d t/3 VC^ 1 + A s ^ 



where 



2 ), 



2a 2a 
A u = -^fi*M;i(M 2 0i ,M 2 ^,\M 3 \ 2 ), A c = ^/i*M*J(M? 2 ,M? 2 ,|M3 

A d = ^/^/(M^M^lMal 2 ), A s = ^ //M* J(M? 2 , M| 2 , |M 3 | 2 ) , (2.20) 



and 



ln(z/y) + yz \n(y/z) + :rz ln(z/a:) 
i[x,y,z) — 

{x -y){y- z) [x - z) 

is the one-loop function that takes account of the non-holomorphic threshold corrections. 



3 Higher-Order Contributions to EDMs 

Beyond the one-loop single-particle level, there are many CP-violating operators that can 
have significant effects on the EDMs. Of particular importance are the gluino- and Higgs- 
mediated dimension-6 Weinberg operator, the Higgs-exchange four-fermion operators, and 
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Figure 2: Typical Feynman diagrams for (d G )Q and (d G ) H . The Hi lines denote all three 
neutral Higgs bosons including CP-violating Higgs-boson mixing. Heavy dots indicate re- 
summation of threshold corrections to the corresponding Yukawa couplings. 

the two-loop Higgs-mediated Barr-Zee-type diagrams. We present detailed analytic ex- 
pressions for all those contributions in this Section. 



3.1 Weinberg Operator 

We start by considering first the gluonic dimension-6 Weinberg operator. This is described 
by the interaction Lagrangian§ [21]: 

^Weinberg = ^ d° f a b c G^ & G° v P , (3.1) 

where G^ = \e^ x<7 G Xa is the dual of the SU(3) C field-strength tensor G Xa . In the MSSM, 
d G is induced at two-loop order by g and t, b, and at three-loop order by CP-violating 
Higgs-boson mixing and non-holomorphic threshold corrections to the top- and bottom- 
quark Yukawa couplings. Hence, the Weinberg operator d G is the sum of two terms: 

d G = (d G )' 9 + (d G ) H , (3.2) 

as illustrated in Fig. [2j The first term, (d G )®, is the quark-squark-gluino exchange contri- 
bution and is given by [19] 

(d G Y 9 = -—(—^—] T m„\y -^Sm[(g^Ygi q f]\ 
V ; 2vr V4vr|M 3 |y q ^ b q \ y |M 3 | 2 Lw-Rj/ ^Lji j 

xH(ml/\M 3 \ 2 ,ml/\M 3 \ 2 ,m 2 q /\M 3 \ 2 ) , (3.3) 

where 

1 f 1 f 1 f 1 NiN 2 

H(z 1 ,z 2 ,z q ) = - dx du dyx(l-x)u—— (3.4) 

2 jo Jo Jo 

*Note the coefficient d G has a different sign from that in [21]: d G = —C. 
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Figure 3: The functional dependence of z g \H(zi, z 2 , z q )\ on z q for several values of z s = z\ 
z 2 . When z s < 4, H(z\, z 2 , z q ) becomes negative when z q is larger than a specific value. 



and 



Ni = u(l — x) + z q x(l — x)(l — u) — 2ux[ziy + z 2 (l 



N 2 
D 



(l_ x )2( 1 _ n) 2 + n 2_^ x 2 (1 _ n) 2 ! 

u(l — x) + Z q x{l — x)(l — U) + Ux[Ziy + Z 2 (l — %))] . 



(3.5) 



Figure [3] shows the functional dependence of z q \H(zi, z 2 , z q )\ on z q for several values of 
z s = Zi = z 2 . Our results are at variance with those presented in [19]. For instance, we 
find that when z s < 4, H(zi, z 2 , z q ) becomes negative beyond certain values of z g §. 



^We thank Oleg Lebedev and Pran Nath for useful comparisons and comments regarding the loop 
function H(z±, z 2 , z q ). 
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The loop function H(zi,Z2,z q ) simplifies considerably in specific regions of the pa- 
rameter space. Specifically, if z\ = z 2 = z s and the limit z q — > is taken, the following 
function may be defined: 



H (z a ) = lim z q H(z„z s ,z q ) 

Zg^O 



18(1 -z s ) 4 



2(1 - z s ){l + Uz a ) - (1 - 16z s - 9z 2 s ) \nz t 



(3.6) 



with H (l) = 5/108. Hence, for z q < 0.1 and (z 2 — 2a) 2 /(£2 + ^i) 2 *C 1, the loop function 
H(zi, z 2, z q) m ay be approximated as follows: 



H(zi, Z2, z q) 



Z\ + z 2 



,Z 2 ~ ZiY — /2i + 2 2 



^1 + z 2 



(3.7) 



where 



108(1 - z* 



1 - z s )(l + 72, + 295^ 2 + 177z 
+6z 2 (21 + 502 s + 9z s 2 )ln(z s ) 



(3J 



with ^(1) = 11/1080. 

The second term in (I3.2p . (d G ) H , is the neutral Higgs contribution [20,21], which may 
be cast into the form: 



(d 



G\H 



4V2G F g 3 s 

(47T) 



q=t,b 



X/ 3^59 fi^gi? 



(3.9) 



where = Mfjjm 2 and 
/i(z) = 



1 r 1 z" 1 m 3 x 3 (1 — x) 

4io Jo [x(l - tii) + z(l - m)(1 - x)} 2 



(3.10) 



with /i(0) = 1/16. We note that for the loop function h(z) we follow [20], whose result is 
smaller by a factor 2 than the one given in [21]. 



3.2 CP-Odd Four-Fermion Interactions 

CP-odd four-fermion interactions play a significant role in the EDMs. These interactions 
may be generically described by the Lagrangian 

Af = XQK/yH/v- :,/'). (3.H) 

/,/' 
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Figure 4: Feynman diagrams for CP-odd four-fermion operators. The Hi line denotes all 
three neutral Higgs bosons including CP-violating Higgs-boson mixing. Heavy dots indicate 
resummation of threshold corrections to the corresponding Yukawa couplings. 

The CP-odd four-fermion operators in ( 13. lip are generated by CP-violating neutral Higgs- 
boson mixing in the t-channel and by CP- violating Yukawa threshold corrections, see Fig. HJ 
The combined effect of these two contributions gives rise to the CP-odd coefficients 



where gj = irif/v with v = 2M\y/g for / = l,d,u. Possible sub-dominant contributions 
from box diagrams [17] have been neglected. 

3.3 Barr— Zee Graphs 

Finally, there are additional Higgs-boson quantum effects that contribute significantly to 
the EDMs beyond the one-loop level. For the Thallium EDM, these are the two-loop 
Barr-Zee graphs, denoted as (df) H , and the CP-odd electron- nucleon interaction Cc s = 
Cs ery 5 e NN [12,13], which is induced by CP-violating gluon-gluon-Higgs couplings, (Cs) 9 - 
As shown in Fig. [51 the electron EDM (d^) H is induced by CP-violating phases of third- 
generation fermions and sfermions and of charginos. More explicitly, (df) H is given by 



(C//0 



H 



9f9f> H 




(3.12) 




1 q 97r 3 2^ n,f2 2^ 9H l q*q ] ? Vqji) 
\ AL1X j=l m Hi j=l,2 



cm ^ q " x - p s ff \ i s p ( \ 

2 a 2 M 2 2^ [dH.e+e-dHiqq J \ T qi) + 9H i e.+e-9H i qq 9{ T qi) 



W 1V1 W i=l 



+ 



+ 




Yl [^e+e-^r+r- f( T ri) + 9 S H i e+e-9H i T+r- 9^) 



10 



>7,0) 



Hi/ 



7 



/ 



rr, (g) 
f 



r,Lb ' 






Figure 5: Barr-Zee diagrams: the Hi lines denote all three neutral Higgs bosons, includ- 
ing CP-violating Higgs-boson mixing, and heavy dots indicate resummation of threshold 
corrections to the corresponding Yukawa couplings. 



+ 



AV2tt 2 s 2 w M w 
3 1 

*EE — 



i=i 3=1,2 m xf 



(3.13) 



with r X i = m 2 x /M'j Ii . The Higgs-mediated two- loop quark EDMs (dq =uds ) H are also cal- 
culated similarly. In the above, the two-loop functions F(t), /(t), and g[r) are given 
by 



F{r) 
f(r) 
9(r) 







f 


dx 


Jo 






f 


2 . 


h) 




f 


2 . 


'0 



x(l — x) 



r — x(l — x) 



In 



x(l — x) 



dx 



2 Jo x(l — x) — r 

dx 



2 Jo x(l — x) — t 
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' x(l — x) 
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x(l — x) 
r 



(3.14) 
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Figure 6: Feynman graphs contributing to (Cs) 9 '■ the Hi lines denote all three neutral Higgs 
bosons, including CP-violating Higgs-boson mixing, and heavy dots indicate resummation 
of threshold corrections to the corresponding Yukawa couplings. 

There are subleading two-loop contributions which we neglect [28,29]. Instead, we consider 
the gluon-gluon-Higgs contribution to Cs, see Fig. [61 This is given by 

(Cs) 9 = (0.1 GeV) ^ £ , (3.15) 

v i=l M H t 

where (N\^G a '^G^ v \N) = -(0.1 GeV)iViV is used, and we use the tree-level value #£ ge = 
— tan/5 O a i. In addition, gfj.g g is the scalar form factor Sf in the heavy (s)quark limit: 

^Efe^-g.E^fl (3-16) 

Q=t,b { 3=1,2 qj ) 

where x t = 1 and x 6 = (1- 0.25k) with « = (iV|m s ss|iV)/220 MeV ~ 0.50 ±0.25 [cf. (IQIV . 

Apart from EDMs, the two-loop Barr-Zee graphs also generate CEDMs for the w and 
d quarks. These additional contributions are given by 

q=t,b \ t=l Hi j = l,2 



9s Cts &cm m qi \ ^ [ p S ft \ , S P / \1 I /n im 

167T 2 S^M^ [9H t q m 9H t qq J\ T qi) + 9 Hm qi 9 H t qq 9 {Tqi ) \ j , l^-J-'J 

with = m, gL 



4 Thallium, Neutron, Mercury and Deuteron EDMs 

In this Section we present analytic expressions for the Thallium, neutron, Mercury and 
deuteron EDMs in terms of the constituent particle (C)EDMs and the coefficients of the 
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dimension-six Weinberg operator and the four-fermion operators. 



4.1 Thallium 

We first consider the atomic EDM d^i of 205 T1. This receives contributions mainly from 
two terms [30,31]: 

d T i[ecm] = -585-df[ecm] - 8.5 x 1(T 19 [ecm] • (C s TeV 2 ) + (4.1) 
where <if is the electron EDM, which is given by the sum 

df = (df)^ + (dff + (df ) H . (4.2) 
The coefficient Cs is calculated as 

C s = {CsY f + (Cs) 9 , (4.3) 
where the (down- type) quark contribution is given by [17] 

m « = C™«» + Cj-^^. . (4.4) 
m d m s 

We neglect the u- and c-quark contributions and absorb the fe-quark contribution into {Cs) 9 
together with the top and heavy-squark contributions [cf. fl3.16p ]. We recall that the ratio 
C q e/m q oc m e /v 2 is scale- invariant [17]. More details concerning the calculations of these 
coefficients are given in Appendix lAl 



4.2 Neutron 

We now turn to the calculation of the neutron EDM, d n . Its magnitude is somewhat 
uncertain, because of non-perturbative dynamics at the hadron level. We consider three 
different hadronic approaches for computing d n : (i) the Chiral Quark Model (CQM), (ii) the 
Parton Quark Model (PQM) and (iii) the QCD sum-rule technique. 

The CQM is a non-relativistic model, where the quark EDMs are estimated via naive 
dimensional analysis (NDA) [32], and the neutron EDM is given by 

dn = - 3 dT A ~ld™\ 

d™ A d = V E df + r, C ^d° + V G ^d G , 

<J3 = (d^ c) r + (< (c) )*° + (V c ¥ + > ( 4 - 5 ) 

where the chiral symmetry breaking scale A ~ 1.19 GeV and the rj E,c,G are QCD correction 
factors that describe the renormalization-group (RG) evolution of d E,c and dP from the 
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electroweak (EW) scale, e.g., the Z-boson mass Mz, down to the hadronic scale. These 
QCD correction factors are given in [33] as 



r] E ~ 1.53 , 7] c ~ r] G ~ 3.4 . (4.6) 

We note that the EDM operators d E,c and d G in (14.51) are computed at the EW scale. 

Another approach to computing the neutron EDM is based on the PQM, which uses 
low-energy data related to the constituent-quark contributions to the proton spin combined 
with isospin symmetry [34]. In the PQM, the neutron EDM is given by 

d n = V E (A^ M df + A^ M d E + A^ M df) , (4.7) 

where we consider the following particular values, to contrast with the PQM predictions 

A^ QM = 0.746 , A^ QM = -0.508 , A^ QM = -0.226 . (4.8) 

As in the CQM, d E 's in (14.71) are evaluated at the EW scale. 

The third approach to computing the neutron EDM employs QCD sum rule tech- 
niques [17,35-38]. With the aid of these techniques, the neutron EDM is determined by 



d n = d n (d E d c ) + d n (d G ) + d n (Q 



bd) 



dn(df,d c ) = (1.4 ± 0.6) (d E - 0.25 d E a ) + (1.1 ±0.5) e(d c + 0.5 d c )/g s 
d n (d G ) ~ ±e(20± 10) MeV d G , 



d n {C bd ) ~ ±e2.6 x 10" 3 GeV 5 



Cbd , n _ e Cdb 
+ 0.75 

m b m b 



(4.9) 



where d E and d c should be evaluated at the EW scale and d G at the 1 GeV scale §. We note 
that the contribution of d G to d n is a factor ~ 2 smaller than in the CQM. We calculate 
the coefficients Cm, an d Cdb by means of ( 13. 12ft and evaluate the coefficients g b and gj, at 
the energy scales nib and 1 GeV, respectively. For defmiteness, in our numerical estimates 
we assume that both d G and Cm contribute positively to d n . 

The comparisons between the results obtained in these three approaches indicates the 
significance of the non-perturbative uncertainties in calculating d n . Related uncertainties 
appear also in the calculations of EDMs of nuclei. In Section [51 we present numerical 
estimates of the neutron EDM based on the CQM, the PQM and QCD sum-rules. 



1 We note that isospin symmetry between the neutron n and the proton p implies that = (A u ) p = 4/3, 
A„ = (A(i)p = —1/3. Furthermore, in the relativistic Naive Quark Model (NQM), one has A s = (A s ) p = 0. 
HHere we make use of the relation: d G \ 1 GoV ~ (?7 G /0.4) d G | EW ~ 8.5<i G | EW [37]. 
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4.3 Mercury 

Using QCD sum rules [17,38] to calculate the Mercury EDM, one finds that 
d Rg = 7x 10~ 3 e{dg-d$)/g a 
- 1.4 x 10" 5 eGeV 2 
+ (3.5 x 1(T 3 GeV)eC, 



2 lE 



h 3.3k h (1 — 0.25k) 

m d m s m b . 



+ (4 x 1(T 4 GeV) e 



s 
C P 



Z -N 
A 



C' 



(4.10) 



The parameters Cp and C' P are the couplings of the CP-odd singlet and triplet electron- 
nucleon interactions, respectively, and are described by the interaction Lagrangian 



C Cp = C P ee Ni^N + C' P ee Ni^ 5 r 3 N . 



(4.11) 



Making use of the SU(2) isospin symmetry of the nucleon N = (p,n), one may evaluate 
the triplet contribution C' P , which is suppressed by a factor [(Z — N)/A] H = —0.2 with 
respect to the singlet one Cp. In detail, the contribution of four-fermion interactions to Cp 
and Cp is given by 



Cr 



C' 



(Cp) 4/ * 



C 



375 MeV J2 ~ 



q=c,s,t,b y 



c 



06 MeV — - 181 MeV £ ^ 



C, 



q=c,s,t,b y 



(4.12) 



where the w-quark contribution has been neglected. For a detailed discussion, see Ap- 
pendix [Bl 

On the basis of above results, we improve upon an earlier calculation of the Mercury 
EDM [310 More explicitly, including the CP-odd triplet electron-nucleon interaction, the 
Mercury EDM is given by 



j nNN 



a {1) 



1.8 x 10~ 3 GeV 



f 10~ 2 df + (3.5 x 10~ 3 GeV) eC s 
Z -N" 



+ (4 x 10" 4 GeV) e 



C P + 



2t\ x 10 



-12 



{d C u ~ d c d )/g s 



A 

(qq)\ 



C' 



Hg 



10- 26 cm (225 MeV) 3 
-8 x 10~ 3 GeV 3 



0-5C dd C sd , 
h 3.3k h (1 

m d m s 



0.25k) 



C 



bd 



(4.13) 



where C nNN D g { ^ N NNn . We note that the factors 1.8 x 10~ 3 GeV" 1 and -8 x 10~ 3 GeV a 
are known only up to 50 % accuracy [39]. 



**We also correct the errors in the contributions of to dn g and of four-fermion interactions to 

& [39]- 
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4.4 Deuteron 



Finally, the deuteron EDM may be calculated by using QCD sum rules. We include the 
contributions from the pion-nucleon-nucleon isospin-triplet coupling g^ N [40], the con- 
stituent proton and neutron EDMs [36], and the dimension-six Weinberg operator [21,37]. 
The deuteron EDM dp is then found to be [6] 



d 



D 



jttNN 
a D 



ittNN 
a D 



d D (d n ,d p ) + d D (d G ) 



-(i) 
eg^NN g^NN 



-(1.3±0.3)es£L 



dD(d n , dp) 
d D (d G ) 



Ylixm^ (1 + 20 2 ~ v— — -y-^iviv 
(0.5 ± 0.3) (df + df ) - (0.6 ± 0.3) e [(d c - d c d )/g s 
d n (d G ) + d p (d G ) ~ ±e(20± 10) MeVd G , 



0.3(d c + d u )/g 



C\ 



(4.14) 



where g^NN — 13.45 and £ = ^/nipe/m^, with e = 2.23 MeV being the deuteron binding 
energy. Collecting all the above intermediate results, we may write dp in the more compact 
form: 



d 



D 



-11 



(0.6 ±0.3) e(d c -d u )/g 



(0.2 ± 0.1) e + + (0.5±0.3)« 



+ (1 ±0.2) x 10- 2 eGeV 2 



d E d ) 

°^ + 3.3k^ + (1- 



0.25k 



Cbd 



± e(20± 10) MeVd 



G 



(4.15) 



In the above, d G is evaluated at the 1 GeV scale, and the coupling coefficients g<i,s,b that 
occur in Cdd,sdM are computed at energies 1 GeV, 1 GeV and m b , respectively. All other 
EDM operators are calculated at the EW scale. We observe that the leading dependence 
of do on d G d is the same as in dn g . In the numerical estimates given in the next Section, 
we assume that d G contributes positively to dp. 



5 EDM Constraints 



In this Section, we present illustrative constraints on key soft SUSY-breaking parameters 
and CP phases in the trimixing, the CPX and the MCPMFV scenarios. We use the following 
current experimental limits on the Thallium [1], neutron [2], and Mercury [3] EDMs: 

|d T i| < 9 x 10~ 25 ecm, 
|d Hg | < 2 x 10~ 28 ecm, 
\d n \ < 3xl0~ 26 ecm. (5.1) 

On the other hand, the projected sensitivity to the deuteron EDM is [4] 

\d D \ < (1 - 3) x 10~ 27 ecm. (5.2) 
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For our numerical study, we take 3 x 10~ 27 e cm as a representative expected value. However, 
we note that much better statistical precision at the level of 10~ 29 ecm may be possible in 
principle [5]. 

5.1 Trimixing Scenario 

The trimixing scenario is characterized by large tan/? and a light charged Higgs boson, 
resulting in a strongly-mixed system of three neutral Higgs bosons with mass differences 
smaller than the decay widths [25]: 

tan/5 = 50, M H ± = 155 GeV, 

Mq 3 = M 03 = M bi = M u = = Msusy = 0.5 TeV, 

|/z|=0.5 TeV, \A t>biT \ = 1 TeV, |M 2 | = |Afi| = 0.3 TeV, |M 3 | = 1 TeV, 

$ p = o , $ 1 = $ 2 = 0°, $ Ar = $ Ae = $ Au = $ Ac = = $ As = 0° . (5.3) 

Note that, in this scenario, only two independent CP-violating phases generate EDMs: 
$A = &A t — &A b arid $3. In addition, we introduce a common hierarchy factor p between 
the masses of the first two and third generations: 

M x^=P M x„ (5-4) 

with X = Q,U, D, L, E. For the size of first two generation A terms, we take \A e \ = \A T \, 
\A U>C \ = \A t \, and \A djS \ = \A b \. Note that in this scenario, only d^ ds and d^ d depend on 
the hierarchy factor p through the resummed threshold corrections (12.191) . 

In Fig. [7J we show the absolute value of the Thallium EDM divided by its current 
experimental limit, in the $3-$a plane (left) and as a function of $3 taking $^ = 60° 
(right). In the left frame, the plane is divided into 4 regions: |g?ti/^ EXP | < 1 (black), 
1 < \d Tl /d EXP \ < 10 (red), 10 < \d T i/d EXP \ < 100 (green), and 100 < \d Tl /d EXP \ (magenta). 
The unshaded region is not allowed theoretically. In the right frame, the constituent con- 
tributions from the electron EDM df and the CP-odd electron-nucleon interaction C$ are 
shown in the thin solid and dashed lines, respectively. The thick solid line is for the total. 
Note the non-trivial cancellation between df and Cs contributions around $3 = 275° in 
the right frame. 

In Fig. [HJ we show the neutron EDM calculated in the CQM for two values of the 
hierarchy factor p: p — 1 (upper) and 3 (lower). In the left frames, the regions are 
shaded as in the Thallium case, Fig. [71 In the right frames, the constituent contributions 
from the quark EDMs d^ d , the quark CEDMs d^ d , the Weinberg operator d G are shown 
in the thin solid, dashed, and dotted lines, respectively. The thick lines are again for 
the total. We note first that d^ d and d^ d do not vanish, even though we are taking 
$1,2 = §a u d = 0° in this scenario. This is because the non-vanishing $3 enters the Yukawa 
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couplings through the threshold corrections and the two-loop Higgs-mediated graphs. The 

EC (~* 

d u d contributions decreases for larger p, whilst the d contribution is independent of p. For 
larger p — 3, the d G contribution is almost dominating. The non-vanishing p-independent 

E C 

dZ'a at $3 = 0°, 180° comes from the two-loop Higgs-mediated diagrams. A non-trivial 
cancellation occurs, for example, between the d c d and d G contributions around $3 = 10° 
in the lower-right frame. 

In Fig. M, we display the neutron EDM calculated in the PQM. The shaded regions 
and lines are the same as in Fig. El except for the constituent contributions in the right 
frames: the thin solid, dashed, and dotted lines are for the contributions from the EDMs of 
u, d, and s quarks, respectively. In this model, the neutron EDM is dominated by the df 
contribution and decreases as p increases. The dips in the lower-right frames are due to the 
cancellation between the contributions to df from the Higgs-mediated and other diagrams. 

Figure fTUl shows the predictions for the neutron EDM calculated using the QCD sum 
rule approach. In the right frames, the thin solid, dashed, dotted, and dash-dotted lines 
are for the constituent contributions df d , d c d , d G , and C^db, respectively. Compared to 
the CQM, d c d contribution is about 3 times larger and d G one about 2 times smaller. The 
Cbd,db contribution is significant because of the large value of tan j3 and the light Higgs 
spectrum in this scenario. 

Comparing Figs. El E] and [101 we see that they yield qualitatively similar overall 
results, but with important detail differences. In particular, the appearances and locations 
of non-trivial cancellations are mo del- dependent. 

Finally, Fig. [11] gives our numerical estimates for the Mercury EDM. In the right 
frames, the thin solid, dashed, dotted, and dash-dotted lines are for the constituent con- 
tributions from df, d c d , C 4 f = C dd>sd ^ d , and C^ P , respectively. When p = 1, the Mercury 
EDM is dominated by d G d . However, as p increases, the d G d contribution decreases while 
other four contributions remain the same. Cancellations occurs more easily for larger p, in 
which case all the contributions become more or less comparable, and we see a non-trivial 
example in the lower-right frame. 

5.2 CPX Scenario 

In the CPX scenario [24], the product of p and the third-generation A terms are larger 
than the common SUSY scale of the third-generation squarks by a factor of 8: 

M Q 3 = M v, = M b, = M i3 = M hi = M SUSY , 

|/i|=4M SU sy, \A tM \ = 2M SUSY , |M 3 | = 1 TeV. (5.5) 

As an example, we have fixed \M 2 \ = 2|Mi| = 100 GeV, and taken the charged Higgs- 
boson pole mass M H ± = 300 GeV and the common SUSY scale M SU sy = 0.5 TeV, but the 
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parameter tan (3 is varied. The A-term phases of first two generations are set to vanish, 
as in the trimixing scenario: $A r = $A e = &A U = &A C = &A d = &A S = 0°. For the 
size of the A terms of the first two generations, we also take \A e \ = \A T \, \A UyC \ = \A t \, 
and \Ad )S \ = \Ab\. As for the CP- violating phases, initially we vary two phases generating 
EDMs, $,4 and $3, taking $1 = $2 = 0. The effects of non-trivial $1 and $2 are described 
later. 

Taking $1 = $ 2 = 0°, in Figs. 112], [TBI [Ti | US ] and[U, we show the Thallium, neutron, 
and Mercury EDMs on the $3 — $a plane (left) and as functions of $3 taking $^4 = 90° 
(right). We take two values of tan/3 = 5 (upper) and 50 (lower) but with fixed p — 1. The 
shaded regions and the lines are the same as in Figs, El EH [TQl and [HI respectively. 

We observe the Thallium EDM, shown in Fig. [12], is dominated by the two-loop Higgs- 
mediated electron EDM for both cases, with only mild dependence on $ 3 . For small tan/3, 
we always have \dTi/d EXP \ < 10, independently of $3,4, as can be seen from the upper-left 
frame. For large tan/3, the sub-leading Cs contribution becomes larger by two orders of 
magnitude, whilst the electron EDM contribution is larger by one order of magnitude as 
seen by comparing the two right frames. 

Turning now to the neutron EDM in the CQM, shown in Fig. [131 we see that the three 
contributions from df d , d^ d , and dF are comparable. The most important contributions 
to the down-quark EDM df and CEDM d c d come from the one-loop gluino diagrams, 
explaining the mild $,4 dependence. The different dependence on $3 for large tan/3 is 

E C 

due to the enhanced two-loop Higgs-mediated contribution to d d ' , and we note a non- 
trivial cancellation in the lower- right panel when $ 3 ~ 325°. The neutron EDM in the 
PQM, Fig. [141 and that found using the QCD sum rules, Fig. [151 are somewhat larger 
in general and show a similar behaviour, due to the dominance by the df and the d d 
contributions, respectively. The Mercury EDM, Fig. [161 is also dominated by the d d 
contribution. However, the sub-dominant contributions from df, C±f = Cdd,sd,bdi anc ^ ^Ip 
become larger for large tan/3. 

We now study the effects of non-zero $1 and $2, and some cancellation properties, 
varying the common hierarchy factor p and assuming maximal CP violation in $^ and 
$3, i.e., <3>a = $3 = 90°. The one-loop contributions decrease as p increases but the con- 
tributions from the Weinberg operator and the two-loop Higgs mediated diagrams remain 
constant. 

In the upper frames of Fig. [T71 we show the Thallium EDM as a function of p for 
two values of tan/3: 5 (upper-left) and 50 (upper-right). The four lines are for ($ 1; $ 2 ) = 
(0°,0°) (solid), (90°, 0°) (dashed), (0°,90°) (dotted), and (0°,270°) (dash-dotted). When 
($1,^2) = (0°,0°), the Thallium EDM is independent of p because the main contribution 
from the electron EDM is dominated by the two-loop Higgs-mediated diagrams, see Fig. [T2l 
The two dips around p = 1.5 and 4 when ($i,$ 2 ) = (90°, 0°) and (0°,90°) are due to 
the cancellations between the p-independent two-loop Higgs and neutralino and chargino 
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contributions to the electron EDM, respectively. In the lower frames we show explicitly the 
chargino-Higgs cancellation in the (0°, 90°) case. In the lower frames, the thick solid line 
is for the total electron EDM and the thin solid, dashed, horizontal dash-dotted lines are 
for the chargino, neutralino, and Higgs contributions to it, respectively. 

In Fig. [TU we compare the three calculations of the neutron EDM: the CQM (upper 
row), the PQM (middle row), and the QCD sum-rule technique (lower row). The lines are 
the same as in Fig. Q3 The cases with $ 2 ) = (0°, 0°) (solid) and (90°, 0°) (dashed) are 
hardly distinguishable from each other. The PQM calculation (middle) is most sensitive 
to $2 ; whilst the QCD sum-rule approach shows the least sensitivity. We observe that the 
value of p where the cancellation occurs varies on the models and approaches. In all cases, 
the neutron EDM calculations saturate to certain values determined by the p-independent 
contribution from the Weinberg operator or the two-loop Higgs mediated diagrams, as 
shown below. 

In Figs. IT9landl2"0l we show details of the neutron EDM taking the cases ($i, $2) = 
(0°,90°) and (0°,270°), respectively, as examples. In Fig. [191 the thick line is for the 
total neutron EDM and the thin solid, dashed and horizontal dotted lines are for the 
contributions from the EDMs of quarks d^ d , the CEDMs d c d , and the Weinberg operator 
d G , respectively. The EDMs and CEDMs are dominated by those of the down quark. In 
the lower-right frame, the lower horizontal dash-dotted line is for the contribution from 
the bottom-down four-fermion operator, Cbd,db- We observe the dips of the thick lines 
are determined by the interplay of the three main contributions. For example, in the 
upper-left frames with tan/3 = 5, cancellation occurs at p = 2 in the CQM. On the other 
hand, in the QCD sum-rule approach where the contribution from d G (d° d ) is suppressed 
(enhanced) compared to the CQM, the cancellation occurs at p ~ 3.6. Taking into account 
of the uncertainty involved in the calculation of the d G contribution, the p value where the 
cancellation occurs may change by ~ ±1, at least. The more significant contribution from 
d^ d in the CQM explains why this calculation is more sensitive to $ 2 than the QCD sum- 
rule approach. The dips in the thin solid lines for d^ d are due to three-way cancellations 
among the one-loop gluino, one-loop chargino and two-loop Higgs-mediated diagrams. For 
large tan /3, no cancellation occurs due to the dominance of d G d . 

Figure [20] displays our numerical estimates for the neutron EDM in the PQM when 
($1, $2) = (0°, 270°). In the upper frames, the thick lines are for the total EDM, and the 
thin solid, dashed, and dotted lines are for the contributions from the EDMs of the up 
d^, down d d , and the strange quark df , respectively. In the lower frame, we show the 
dominating contribution from the strange-quark EDM df (thick solid) together with its 
constituent contributions from the one- loop chargino x (thin solid), one-loop neutralino x° 
(thin dashed), one-loop gluino g (thin dotted), and two-loop Higgs-mediated H° (horizontal 
thin dashed-dotted) diagrams. We observe that there is a cancellation between the chargino, 
gluino, and Higgs contributions at p ~ 4. 
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The above examples demonstrate that the interpretation of neutron EDM measure- 
ments is subject to uncertainties in non-perturbative QCD, which are reflected in the 
specific models discussed. 

In Fig. [2TJ we show the Mercury EDM. It shows barely any sensitivity to $1 and $2, 
due to the dominance by the down-quark CEDM d c , which is dominated by the one-loop 
gluino and two-loop Higgs-mediated diagrams, as shown in the lower frames. In the lower 
frames, the lines are the same as in the lower frames of Fig. 1201 

Taking into account the uncertainty of the d G contribution to the neutron EDM, 
and that involved in the Mercury EDM calculation, there is a possibility of evading all 
the three EDM constraints by taking ($1,^2) ~ (0°,90°) and p ~ 4 in the CPX scenario 
with $^4 = $3 = 90° when tan (3 = 5, if one calculates the neutron EDM using the QCD 
sum-rule approach. 

Finally, in Fig. [22] we show the deuteron EDM. It is not sensitive to $1 and $2; as 
it is dominated by the d G and d c contributions. At high p, the EDM is given by the 
sum of the dF contribution and the contribution to d c from the two-loop Higgs-mediated 
diagrams. In the lower frames, the thick line is for the total EDM and the thin solid, 
dashed, horizontal dashed, horizontal dash-dotted lines are for the contributions from d^ d , 



5.3 MCPMFV Scenario 



In the MCPMFV scenario, in contrast to the trimixing and CPX scenarios, the soft SUSY- 
breaking parameters are specified at the the gauge coupling unification (GUT) scale where 
the MFV condition is imposed [26,41]. This scenario has a total of 19 parameters, including 
6 CP-violating phases and 13 real mass parameters. As a numerical example, in order to 
study the effects of the CP-violating phases on EDMs in this framework, we consider a 
CP-violating variant of a SPSla-like [42] scenario: 

|Afi l2>3 | = 250 GeV, 

^=^ = ^ = ^ = ^ = (100 GeV) 2 , 
M 2 = M| = (200 GeV) 2 , 

\A U \ = \ A d \ = \A e \ — 100 GeV, (5.6) 

with tan/3 (Msusy) = 40. As for the CP- violating phases, we adopt the convention that 
$ M = 0°, and we vary separately the three phases of the gaugino mass parameters, $1,2,3, 
taking however vanishing A-term phases at the GUT scale: = &a^ T = ^ < A^ T = 0°- 

At the low-energy Msusy scale, A-term phases may be generated by the CP- violating phases 
of the gaugino mass parameters $1,2,3, even though we are taking real A terms at the GUT 
scale [26,43]. 
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In the upper-left frame of Fig. [23j we show the Thallium EDM as a function of $2 
for several values of ($i,$ 3 ): (0°,0°) (solid), (0°, 90°) (dashed), (90°, 0°) (dotted), and 
(270°, 0°) (dash-dotted). We see the Thallium EDM is nearly independent of $3, cf., the 
two overlapping lines for ($ 1; $ 3 ) = (0°,0°) and (0°,90°). This is because the Thallium 
EDM is dominated by the electron EDM. In the upper-right, lower- left, and lower-right 
frames, we show the electron EDM as a function of $ 2 when ($1, <3> 3 ) = (0°, 90°), (90°, 0°), 
and (270°, 90°), respectively. In the lower frames, we observe the neutralino contribution 
becomes less dependent on $ 2 when $x = 90° and 270° and cancels the chargino one around 
$2 = ±4°, resulting in the dips in the upper-left frame. 

Figure [2H presents numerical estimates of the neutron EDM in the CQM as a function 
of $2- We observe that the EDM is nearly independent of $1 and $3. This is because of 
the dominance of the down-quark EDM d d , see the upper-right frame with ($1,^3) = 
(0°,0°) in which the thin solid line for d^ d is overlapped by the total thick line. The 
electric EDM of the down quark still plays a dominant role even when $ 3 has a non-trivial 
value, because of an accidental cancellation between the contributions from the Weinberg 
operator d G (dotted line) and the down-quark CEDM d d (dashed line) as shown in the lower 
frames for ($i,$ 3 ) = (180°, 90°) (lower-left) and (180°, 270°) (lower-right). This accidental 
cancellation is lifted when the EDM is calculated using the QCD sum-rule approach, in 
which d G and d c contribute differently to the total EDM, see the lower-frames of Fig. [261 

In the upper-left frame of Fig. [251 we show the neutron EDM calculated in the PQM. 
We find that the neutron EDM is nearly independent of $ 3 . The dominant contribution 
comes from the strange-quark EDM df , and we show it together with the constituent 
contributions in the other frames taking ($i,$3) = (0°,0°) (upper-right), (0°, 10°) (lower- 
left), and (0°,90°) (lower-right). The thick lines are for the total EDM and the thin solid, 
dashed, dotted, and dash-dotted lines are for the contributions from the one-loop chargino 
X ± , neutralino \°i gluino g, and the two-loop Higgs-mediated H° diagrams. When $3 = 0°, 
df is dominated by the chargino contribution, see the upper-right frame. As shown in 
the lower frames, the $3 dependence of the strange-quark EDM is largely cancelled in 
the two main chargino and gluino contributions, explaining the $3 independence of the 
neutron EDM shown in the upper-left frame. We note that this cancellation resulting in 
$3 independence also occurs in the down-quark EDM d d , see the thin solid d^ d lines in 
the lower frames of Figs. [M] and [261 

Figure [26] gives the predicted values for the neutron EDM calculated using the QCD 
sum- rule approach. In the upper-right and lower frames, the lines are the same as in Fig. [2U 
but for different combinations of $1 and $3. The non-trivial three-way cancellation among 
the contributions from d^ d , d c d , and d G explains the dips in the upper-left frame when 
$ 3 = 20° and 340°. 

In Fig. [271 we present numerical estimates for the Mercury EDM. In the upper-right 
and lower frames, we see that the dominant contributions come from the electron EDM 
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df and the down-quark CEDM . The cancellation between them results in the dips in 
the upper-left frame. Note that $3 = 90° could generate an EDM larger than the current 
experimental limit by a factor ~ 400, nearly independent of $2; see the nearly horizontal 
thin solid line in the upper-left frame. 

Finally, Fig. [28] shows theoretical predictions for the deuteron EDM. The deuteron 
EDM is dominated by the down-quark CEDM d^. The sub- leading contributions are from 
the Weinberg operator d G and the down-quark EDM (if. In the upper-right and lower 
frames, we show the constituent contributions to d^. The chargino and neutralino contri- 
butions are dominant when $ 3 = 0, see the upper-right and lower-left frames. However, 
as $3 grows, the gluino contribution rapidly increases, see the lower-right frame. We note 
that the dip around $2 = —20° in the upper-left frame when ($1, <J> 3 ) = (0°, 10°) is due to 
a cancellation between the d^ and <if contributions. 

6 Conclusions 

We have performed a fully-hedged analysis of the Thallium, neutron, Mercury and deuteron 
EDMs within the general CP-violating framework of the MSSM. In our analysis, we have 
taken into account the the complete set of one-loop graphs, the dominant Higgs-mediated 
Barr-Zee diagrams, the complete CP-odd dimension-six Weinberg operator and the Higgs- 
mediated four-fermion operators. Our study has also improved earlier calculations in two 
important aspects. First, it includes CP- violating Higgs-boson mixing effects and, secondly, 
it properly implements resummation effects due to threshold corrections to the Yukawa 
couplings of all up- and down-type quarks and charged leptons. Not only do these two 
effects turn out to be significant for the existing one- and two-loop EDMs, but also they 
give rise to additional higher-order contributions within the MSSM, such as the original 
Weinberg operator induced by t and b quarks with Higgs bosons in the loop. In addition, 
we have improved the Mercury EDM calculation by including the contribution due to the 
CP-odd triplet electron-nucleon interaction. 

Having established the latest state-of-art theoretical framework as described above, we 
have then explored the EDM constraints on the CPX, trimixing and MCPMFV scenarios. 
Clearly, sufficiently small values of the CP-violating parameters survive the experimental 
upper limits on the EDMs. However, we have also found that larger values of the CP- 
violating parameters may be allowed exceptionally by accidental cancellations among the 
CP- violating contributions to the three measured EDMs, i.e., c£ti, d n and c?Hg, in all the 
above scenarios. In detail, the results of our analysis may be summarized as follows. 

First, we studied the trimixing scenario. We have explored the impact on the three 
measured EDMs resulting from an hierarchy of the soft SUSY-breaking masses between the 
first two and third generations, p, and the CP phases <$>a = = $A b and $3 related to 
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the third generation A-terms and the gluino mass parameter, respectively. In the case of 
the Thallium EDM, large CP-violating phases are allowed due to the cancellation between 
the contributions from the CP-odd electron-nucleon interaction Cs and the two- loop Higgs- 
mediated electron EDM (df) H . The two contributions are found to be independent of p in 
the scenario under consideration. 

The neutron EDM has been computed within three different models and approaches: 
(i) the CQM, (ii) the PQM and (iii) a QCD sum-rule approach. In the CQM and QCD 
sum-rule approaches, the neutron EDM is dominated by the EDM and CEDM of the down 

EC c 1 

quark d d ' and the contribution from the dimension-six Weinberg operator d . In the QCD 
sum rule approach, d G (d d ) becomes less (more) important than in the CQM. The down- 

E C 

quark EDM and CEDM d d ' are generated through the resummed threshold corrections 
at one loop by $3 and the two-loop Higgs-mediated diagrams. The latter remains the 
same for large p. The df' c contributions depend on p, whilst the one due to d° is p- 
independent. The bottom-down four-quark interaction Cu,db is a ls° independent of p and 
becomes significant in the QCD sum rule approach. Interestingly, cancellations may occur 

EC c 

among d d ' , d and Cbd,db more easily when p is large. In the PQM, on the other hand, the 
neutron EDM is dominated by the strange-quark EDM df, which is generated by one- loop 
threshold corrections at the one-loop level and by the two-loop Barr-Zee graphs. Again, 
cancellations occur between the two contributions when p is large. All the three different 
hadronic models or approaches for computing the neutron EDM give comparable estimates. 
However, there are significant differences of detail, implying that the interpretation of 
neutron EDM measurements is mo del- dependent. 

As far as the Mercury EDM in the trimixing scenario is concerned, the <i-quark CEDM 
d d was found to be dominant for small values of p, whilst it gets suppressed for large values 
of p. In addition, the p-independent contributions from Cs, the four-quark operators 
Cdd,sd,bdi the CP-odd singlet and triplet electron-nucleon interactions, Cp and C' P , and the 
electron EDM become all relevant, leading to an interesting cancellation pattern. 

In the CPX scenario, we have analyzed the dependence of EDMs on the hierarchy 
factor p, the CP phases $a = <&A t = &A b and $1,2.3, for a relatively low and large value 
of tan/3: tan/3 = 5 and 50. In the absence of any mass hierarchy between the first two 
and third generation of squarks and sleptons, i.e., p — 1, and for $! = $ 2 = 0°, we 
have found that the Thallium EDM is dominated by the electron EDM, whilst Cs starts 
becoming important only for tan/3 S 10. Correspondingly, the neutron EDM receives 

EC C 

the largest contribution from d d ' and d in the CQM and the QCD sum-rule approach. 
The main contributions to d d ' c come from the one-loop gluino diagrams and the two-loop 
Higgs-mediated diagrams. The latter become important when tan/3 is large. In the PQM, 
the neutron EDM mainly results from df. Finally, the Mercury EDM receives its biggest 
contribution from d d , whilst the contribution due to df becomes important only for large 
tan (3. 
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Varying p and the CP phases $1 and $2, we have found that several cancellations can 
occur within the CPX scenario. Specifically, one-loop neutralino and chargino effects and 
the two-loop Barr-Zee graphs may add up destructively and suppress the electron EDM 
df. Likewise, cancellations among the one-loop chargino and gluino graphs and two-loop 
Higgs-mediated diagrams lead to suppressed strange- and down-quark EDMs df s and to 
an equally small down-quark CEDM d c . As a consequence of the suppressed df, df and 
d° , the Thallium, neutron (in the PQM) and Mercury EDMs were all found to come out 
well below their experimental limits. A similar result has been obtained in the CQM or 
QCD sum-rule approach, where the Weinberg operator d G plays an important role. In 
particular, we have demonstrated explicitly the possibility of evading all the three EDM 
constraints from g?ti, d n and dn g in the CPX scenario with <&a ~ $3 ~ 90°, if one assumes 
that ($i,5>2) ~ (0°,90°) and p ~ 4 when tan/3 = 5, provided d n is calculated using QCD 
sum- rule techniques. Finally, the deuteron EDM [4,5], whose size crucially depends on dP 
and d c , can constrain dramatically the CP- violating phases by a factor of about 100 [6]. 

A third benchmark scenario of the MSSM that has been analyzed was the MCPMFV 
scenario. For definiteness, we have used SPSla-like input parameters given at the GUT 
scale, for which $9 UT = $9 UT = $9 UT = 0° and tan/3 = 40. This choice results in a 
relatively light SUSY spectrum, where phases for the A-terms are generated by RG running 
down to the electroweak scale. Within this particular MCPMFV scenario, we have studied 
the implications of the CP- violating phases $1, $2 and $3 for the three EDMs d^i, d n and 
c?Hg- The Thallium EDM was found to be nearly independent of $3, since the electron 
EDM df turns out to be the dominant contribution. Moreover, we have noticed that one- 
loop chargino and neutralino effects may cancel each another in df. The neutron EDM 
in the CQM was also found to be nearly independent of $1 and $3, due to an accidental 
cancellation between the Weinberg operator d G and the down-quark CEDM d c . Unlike 
in the CQM, this cancellation is no longer present when the QCD sum rule approach is 
used, leading to a strong dependence on $ 3 . In the PQM, like in the CQM, the neutron 
EDM does not depend strongly on $3, since there is a cancellation between the chargino 
and gluino effects on the dominant strange-quark EDM df. The Mercury EDM dn g in the 
MCPMFV scenario was found to receive its biggest contribution from the electron EDM 
df and the down-quark CEDM d c . A non- vanishing gluino phase $3 can drive dn g to quite 
large values that could easily exceed the current experimental limit. Finally, the leading 
effect on the deuteron EDM comes from d c , whilst d G and df remain sub-leading. Since 
the deuteron EDM is expected to be a factor ~ 100 more sensitive than dn g [4,5], it will 
lead to much tighter constraints on the CP- violating phases [6]. 

Our detailed study has shown that the three measured EDMs provide correlated con- 
straints on the 6 CP-violating phases in MCPMFV scenario, leaving open the possibility of 
relatively large contributions to other CP- violating observables. In particular, the deuteron 
EDM [4, 5] will probe the unconstrained CP phases of the MCPMFV scenario. 
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The analytic expressions for the EDMs are implemented in an updated version of 
the code CPsuperH2.0. This new feature of CPsuperH2.0 will be particularly valuable in 
future explorations for possible new-physics phenomena in the K- and 5-meson systems. 
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A Calculation of Cs 



In this appendix we calculate the coefficient C$- Our starting point is the interaction 
Lagrangians: 



C-Cs = C s eij 5 e NN , £ 4f = C qe qq eij 5 e . 
Given the relation (N\£ Cs \N) = (N\£tf \N), one may identify 



c 



qe 



(N\qq\N) 
(N\NN\N) 



(a.i; 



(A.2) 



where q could be a light quark, e.g. u, d, or a heavy one such as the b quark. 



We first derive the light-quark contribution to (Cs) 4 ^ , (Cs) 4 ^ 
need to know (N\uu\N) and (N\dd\N). Using the relation [17] 

(m u + m d ) (N\uu + dd\N) ~ 90 MeV (N\NN\N) 

and assuming that the triplet contribution vanishes, i.e. 

(N\uu - dd\N) = 0, 

we obtain 



q=u,d 



_. To this end, we 
(A.3) 

(A.4) 



(N\uu\N) 



(N\dd\N) 



1 90 MeV 



29 MeV 



(N\NN\N) (N\NN\N) 2 (m u + m d ) m d 
with m u /m d = 0.55. Putting everything together, we find that 



m d 



29 MeV 



(A.5) 



(Cs 



m d 



29 MeV 



16 MeV 



m,. 



(C S ) 4f 



a 



29 MeV 



tie 



(A.6) 



d m d 

The very last expression was used to obtain the first term in Eq. (14.41) . 

As for the heavy-quark contribution to Cs, there are two approaches that can be 
considered. We illustrate these by taking the 6-quark as an example. The first way is to 
include the heavy-quark contribution directly to (Cg) 4 -^ [17], i.e. 

66MeV(l - 0.25k) 



(Cs) 4/ 



a 



be' 



(A.7) 



6 rrib 

The second method uses the QCD trace anomaly and the heavy quark is integrated out in 



the gluon-gluon-Higgs vertex (Cs) 9 [cf. A3. 15j) 

[9h, 



l(Cs) 9 } b 



(O.lGeV)^E 



.</.'/ 



9Hiee 



i=l 



Ml 



a 



be 



2( 0.1GeV) 
3m b 



(A.8) 



where we made use of the relation [g s H . gg ]b = 2/3 g^. u [cf. (l3~T6l) and (l3~T2l) ]. Notice that, 
apart from the ^-dependent term, the two approaches are equivalent to each other. 
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B Calculation of Cp and C' } 



Here we compute the iso-scalar and iso-triplet coefficients Cp and C' P that are relevant in 
the determination of the Mercury EDM. To this end, we start considering the interaction 
Lagrangians, 

£c P = CpeeNij 5 N + C P eeNi^f 5 T 3 N , £ 4 f = C eq eeqi~f 5 q . (B.l) 

Imposing the relation (N\Cc P \N) = (N\C4f\N), we may project out the iso-scalar and 
iso-triplet contributions as follows: 



(N\qij B q\N) r,„, . 4fl {N\qi^ 5 q\N) 



l ° pj J, ~ ° e9 (iV|M 75 iV|iV) ' H CpJ \ q - Ueq (N\m l5 T 3 N\N) ■ { ^ Z) 

As in Appendix A, we need to consider the light- and heavy-quark contributions 
separately. Our approach closely follows [44]. Thus, taking into account all the relations 
that follow from isospin invariance, i.e. 

{N\m u ui l5 u + m d dij 5 d\N) = m N (-g A ) (N\Ni l5 r 3 N\N) , 

m u + m d 

(N\uij 5 u\N) = -(N\dij 5 d\N) , (B.3) 

we obtain 

Tfl — 

m u (N\m lb u\N) = m N (-g A ) f— (N\Ni lh r 3 N\N) , 

m u + m d 

m d (N\di l5 d\N) = -m N (-g A ) ™ d (N\Ni l5 r 3 N\N) , (B.4) 

m u + m d 

where (—gA) — 1-25 is the axial nucleon form factor. From all the above relations, it is 
then not difficult to derive that 

/^ M 4/l c eu , v m u 444 MeV 

(Cp) J = m N (-g A ) ■ ~ C eu , 

u m u m u + m d m u 

<ni \4/l C ed i x rn d 806 MeV 

(C P 4/ = m N (-g A ) : ^ -C ed , B.5 

m d m u + m d m d 

where we assume that = 1 GeV and m u /m d = 0.55. Observe that the light quarks do 
not contribute to the singlet coefficient (Cp) 4 A 

To calculate the heavy-quark contributions to Cp and C' P , we first consider the chiral 
anomaly relations [44]: 



(N^J^N) = 2(N\m u uvy 5 u + m d dij 5 d\N) 

a 
Ttt 



+2 J2 (N\m q qi l5 q\N) + 6(N\^GG\N) 

q=c,s,t,b 



(N\d»J*\N) = (-gf ] )2m N (N\m l5 N\N), (B.6) 
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where rii = 2 and rih = 4 are assumed to be the numbers of the light and heavy 
respectively, and = (3/5)gA in the relativistic quark model. After integrating 
heavy quarks by employing the relation 



(N\m q qi lb q\N)\ 



we get 



q=c,s,t,b 



I n 

—AN\^GG\N) 

A 07T 



(N\m u ui l5 u + m d di l5 d\N) + (N\^GG\N) = i-gf) m N (N\Ni l5 N\N) . (B.8 



quarks, 
out the 



(B.7) 



57T 



From (IB. 31) and (1B.7j) . we finally obtain 



(N\m q qi l5 q\N) = -\m N (-gf) (N\Ni l5 N\N) 



1 m u m d ^ {N \m l5 r 3 N\N) 

2 m u + m d 

for each of the heavy quarks q — c,s, t, b. The first and second terms in (1B.9|) 
heavy-quark contributions to (Cp) 4 -^ and (C P ) 4 ^, respectively. More explicitly, we 



(C P ) 4/ 



(Cp) 4/ 



a 



q=c,s,t,l 



q=c,s,t,l 



eg 



- 2 m N 



(o) 
-9 a 



■a 



375 MeV 



eq 



a 



eq 



1 m u - m d 

2 m u + m d 



m N {-9a) 



-G 



181 MeV 



eq 



(B.9) 

give the 
have 



(B.10) 



where = 1 GeV and m u /m d = 0.55 were used in our numerical estimates. 



C CPsuperH2 . Interface 

• Input: For the complex A parameters of first two generations, part of auxiliary array 
CAUX_H is used as 

A e = CAUX_H(995), 

A u = CAUX_H(996) , A c = CAUX_H(997) , 

A d = CAUX_H(998) , A s = CAUX_H(999) . (C.l) 

• Output: For output, part of auxiliary array RAUX_H is used. 

— The electron EDM in units of cm: 

RAUX_H(200) = df/e = (df /e)** + (df /ef° + (df /ef + (df Je) H , (C.2) 

where the constituent contributions are 

RAUX_H(201) = (df/e)^ , RAUX_H(202) = (df/e)*° 
RAUX_H(203) = (df/e)' 9 , RAUX_H(204) = (df /e) H . 
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— The electric EDM of the u quark in units of cm: 

RAUX_H(210) = df/e = (df/ef* + (d E Jef° + (df/e) 9 + (df/e) H , (C.3) 

where the constituent contributions are 

RAUX_H(211) = (df/e)** , RAUX_H(212) = (df/e)*° 
RAUX_H(213) = (df/e) 9 , RAUX_H(214) = (df /e) H . (C.4) 

— The electric EDM of the d quark in units of cm: 

RAUX_H(220) = df/e = (df/e)** + (df /e)*° + (df/e) 9 + (df/e)" , (C.5) 

where the constituent contributions are 

RAUX_H(221) = (df/e)** , RAUX_H(222) = (df/e)*° 
RAUX_H(223) = (df/e) 9 , RAUX_H(224) = (df /e) H . (C.6) 

— The electric EDM of the s quark in units of cm: 

RAUX_H(230) = df /e = (df /e)** + (df /e)*° + (df/e) 9 + (df /e) H , (C.7) 

where the constituent contributions are 

RAUX_H(231) = (df/e)** , RAUX_H(232) = (df/e)* 
RAUX_H(233) = (df/e) 9 , RAUX_H(234) = (df /e) H . (C.8) 

— The chromo-electric EDM of the u quark in units of cm: 

RAUX_H(240) = d C u = (d° u )** + (d c u )*° + (d C u f + (d C u ) H , (C.9) 

where the constituent contributions are 

RAUX_H(241) = (d%)** , RAUX_H(242) = (d^)*° 

RAUX_H(243) = (d^) 9 , RAUX_H(244) = (d%) H . (CIO) 

— The chromo-electric EDM of the d quark in units of cm: 

RAUX_H(250) = d C d = (d C d )** + (d c d )*° + (d c d ) 9 + (d C d ) H , (C.ll) 

where the constituent contributions are 

RAUX_H(251) = (d%)** , RAUX_H(252) = (d%)*° 

RAUX_H(253) = (d c d ) 9 , RAUX_H(254) = (d c d ) H . (C.12) 
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— The purely gluonic dimension-six Weinberg operator in units of cm/GeV: 

RAUX_H(260) = d G = {d G ) H + (d G Y , (C.13) 

where the constituent contributions are 

RAUX_H(261) = (d G ) H , RAUX_H(262) = {d G )~ 9 . (C.14) 

In the distributed version, Eq. (13. 7p is used to evaluate the function H(zi, Z2, z q ) 
for (d G y . For a full calculation, especially when z q > 0.1, the user should provide 
a dedicated routine, see Fig. [3J 

— Cs, Cp, and, C' P in units of cm/GeV: 

RAUX_H(270) = C s , RAUX_H(27l) = C P , RAUX_H(272) = C' P . (C.15) 

— The coefficients of four-fermion operators in units of cm/GeV 2 : 



RAUX_H(280) 


= C de / m d , 


RAUX_H(281) 


— C se /m s , 




RAUX_H(282) 


= C ed /m d , 


RAUX_H(283) 


= C es /m s , 


RAUX_H(284) = C eb /m b , 


RAUX_H(285) 


= C ec /m c , 


RAUX_H(286) 


= C et /m t , 




RAUX_H(287) 


= Cm/ ™d , 


RAUX_H(288) 


= C sd / m s , 




RAUX_H(289) 


= C bd /m b , 


RAUX_H(290) 


= C db /m b . 


(C.16) 



— The Thallium EDM in units of e cm: 

RAUX_H(300) = cZ T i = d T i(df) + d Tl {C s ) (G.17) 
where the constituent contributions are 

RAUX_H(301) = dn(df ) , RAUX_H(302) = d T i(C s ) . 

— The neutron EDM in units of e cm: 

* Chiral quark model: 

RAUX_H(310) = d n = d n (d^ d ) + d n (d G d ) + d n (d G ) (C.18) 
where the constituent contributions are 

RAUX_H(311) = d n (d^ d ), RAUX_H(312) = d n (d C d ) . RAUX_H(313) = d n (d G ) . 

* Parton quark model: 

RAUX_H(320) =d n = d n (d^) + d n {d^) + d n {df) (C.19) 
where the constituent contributions are 

RAUX_H(321) = d n (d%) , RAUX_H(322) = d n (d d ) . RAUX_H(323) = d n (df) . 
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* QCD sum rule approach: 

RAUX_H(330) = d n = 4(< d ) + d n (d G d ) + d n (d G ) + d n {C bd4b ) (C.20) 

where the constituent contributions are 

RAUX_H(331) = d n (d^ d ) , RAUX_H(332) = d n (d G d ) , 
RAUX_H(332) = d n (d G ) , RAUX_H(334) = d n (C bd ,bd) ■ 

— The Mercury EDM in units of e cm: 

RAUX_H(340) = d Hg = d Hg (df) + d Hg (d c d ) + d Hg (C 4f ) + d Hg (C s ) + d^C^) 

(C.21) 

where the constituent contributions are 

RAUX_H(341) = d ng (df ) , RAUX_H(342) = d Rg (d G d ) , 

RAUX_H(343) = d H g(C 4 /) , RAUX_H(344) = d Kg (C s ) ■ RAUX_H(345) = rf Hg (C?) . 

— The deuteron EDM in units of e cm: 

RAUX_H(350) =d D = d D (d^ d ) + d D (d%) + d n (C 4f ) + d n (d G ) (C.22) 

where the constituent contributions are 

RAUX_H(351) = d D (d^ d ) , RAUX_H(352) = d D (d c d ) , 
RAUX_H(353) = d D (C if ) , RAUX_H(354) = d D (d G ) . 

• IFLAG_H(18) = 1 is used to print out the Thallium, neutron, Mercury and deuteron 
EDMs. Using the run shell-script file distributed, the sample output obtained with 
$A e = $a t , ®a u = ®a c = ®A t , and <$> Ad = $ Aa = $ Ab , and the hierarchy factor 
Pq = Pu = Pd = PL = Pe = 1 is 



Thallium EDM in units of [e cm]: d~Tl/[e cm] 



d~Tl/(e cm) [Total] = -.9371E-24 
Each contribution to d~Tl from 
[d"E_e] = -.9082E-24 
[C_S ]= -.2890E-25 



Neutron EDM in units of [e cm]: d~n/[e cm] 
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(1) Chiral Quark Model 

d~n/(e cm) [Total] = -.4196E-23 
Each contribution to d~n from 

[d~E_u & d~E_d]= -.3994E-23 
[d~C_u & d~C_d] = -.2066E-23 
[ Weinberg-6D ]= 0.1863E-23 

(2) Parton Quark Model 

d"n/(e cm) [Total] = 0.1278E-22 
Each contribution to d~n from 

[d~E_u ]= -.6868E-25 

[d~E_d ]= -.2209E-23 

[d~E_s ]= 0.1506E-22 

(3) QCD sum rule technique 
d"n/(e cm) [Total] = -.7240E-23 
Each contribution to d~n from 

[d~E_u & d~E_d] = -.2741E-23 
[d~C_u & d"C_d]= -.5483E-23 
[ Weinberg-6D ]= 0.9836E-24 
[ C_bd & C_db ]= 0.1749E-28 



Mercury EDM in units of [e cm]: d~Hg/[e cm] 



d~Hg/(e cm) [Total] = 0.3383E-25 
Each contribution to d~Hg from 

[d"E_e ]= 0.1553E-28 

[d~C_u & d"C_d]= 0.3381E-25 
[C_4f ]= -.1094E-29 

[C_S ]= 0.2348E-29 

[C_P & C_P"pr ]= 0.2229E-29 



Deuteron EDM in units of [e cm]: d"D/[e cm] 



d~D/(e cm) [Total] = -.2598E-22 
Each contribution to d~D from 

[d"E_u & d~E_d]= -.9236E-24 
[d"C_u & d~C_d]= -.2604E-22 
[C_4f ]= 0.7811E-27 

[ Weinberg-6D ]= 0.9836E-24 
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IFLAG_H(18)=2 is used to print out the EDMs of the electron and the up, down, and 
strange quarks, the CEDMs of the up, down, and strange quarks, etc. 



The Electric EDMs of particles in cm: e, u, d, s: 



d~E_e/e [Total] 
d~E_u/e [Total] 
d~E_d/e [Total] 
d~E_s/e [Total] 
d~E_e/e[C,N,Gl,H] 
d-E_u/e[C,N,Gl,H] 
d~E_d/e[C,N,Gl,H] 
d~E_s/e[C,N,Gl,H] 



0. 1553E-26 
0.8836E-25 
-. 1936E-23 
-.4355E-22 
0.0000E+00 
0.4467E-28 
0.9825E-25 
0.2211E-23 



-.2833E-26 
0.1755E-26 
-.3481E-26 
-.7832E-25 



0.0000E+00 
0.8717E-25 
-.2042E-23 
-.4595E-22 



0.4386E-26 
-.6115E-27 
0. 1170E-25 
0.2631E-24 



The Chromo-Electric EDMs of particles in cm: u, d: 



d~C_u [Total] 
d~C_d [Total] 
d~C_u [C.N.Gl.H] 
d~C_d [C.N.Gl.H] 



-.1208E-24 
-.5756E-23 

-.8956E-28 0.3073E-26 

0.1626E-25 0.1218E-25 



-.1180E-24 -.5825E-26 
-.5554E-23 - . 2306E-24 



Purely-gluonic D-6 Weinberg operator in cm/GeV: 



d"G [Total]: 0.5786E-23 

d~G[Higgs,Gluino] : 0.2046E-26 0.5784E-23 



Four-fermion couplings needed for EDMs in cm/GeV~2: 



C4_de/m_d 
C4_se/m_s 
C4_ed/m_d 
C4_es/m_s 
C4_eb/m_b 
C4_ec/m_c 
C4_et/m_t 
C4_dd/m_d 
C4_sd/m_s 



-.2432E-26 
-.2432E-26 
-.2432E-26 
-.2432E-26 
-. 1547E-25 
-.9727E-28 
0.3965E-27 
-.3565E-25 
-.3565E-25 
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C4_bd/m_b: 0.1769E-24 
C4_db/m_b : - . 2269E-24 



C_S, CLP, and C_P~prime in cm/GeV and in l/GeV~2: 



C_S 
C_P 

C_P~prime 



0.6710E-27 
0.6602E-26 
0.5147E-26 



0.3400E-13 
0.3346E-12 
0.2608E-12 
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Figure 7: The absolute value of the Thallium EDM divided by its current experimetal limit, 
<if xp = 9 x 10~ 25 e cm, in the $3-$^ plane (left) and as a functon of $3 taking $^ = 60° 
(right). The trimixing scenario has been taken. In the left frame, the plane is divided into 4 
regions: \d T \/d EXP \ < 1 (black), 1 < |d T i/d EXP | < 10 (red), 10 < \d Tl /d EXP \ < 100 (green), 
and 100 < |c?ti/^ exp | (magenta). The unshaded region is not allowed theoretically. In the 
right frame, the constituent contributions from c?f and C$ are shown as the thin solid and 
dashed lines, respectively. The thick solid line is for the total EDM. 
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Figure 8: The absolute value of the neutron EDM in the CQM divided by its current 
experimetal limit, d^ xp = 3 x 10~ 26 e cm, in the $3-$a plane (left) and as a function o/$ 3 
taking $^ = 60° (right). The trimixing scenario has been taken with the common hierachy 
factor p = 1 (upper) and 3 (lower). In the left frames, the shaded regions are the same 
as in Fig. [7| In the right frames, the constituent contributions from d^ d , d^ d , and d G are 
shown in the thin solid, dashed, dotted lines, respectively. The thick solid line is for the 
total EDM. 
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Figure 9: The same as in Fig. [21 but using the PQM for the calculation. In the right 
frames, the constituent contributions from d®, d®, and df are shown as the thin solid, 
dashed, dotted lines, respectively. The thick solid line is for the total EDM. 
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Figure 10: The same as in Fig. [21 but using the QCD sum rule approach for the calculation. 
In the right frames, the constituent contributions from d^ d , d^ d , d G , and Cbd,db ore shown 
as the thin solid, dashed, dotted, dash-dotted lines, respectively. The thick solid line is for 
the total EDM. 
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Figure 11: The absolute value of the Mercury EDM divided by its current experimetal limit, 
^Hg P = 2 x 10~ 28 e cm, on the $3-$ a plane (left) and as a function of $3 taking $^ = 60° 
(right). The trimixing scenario has been taken with the common hierachy factor p — 1 
(upper) and 3 (lower). In the left frames, the shaded regions are the same as in Fig. [7| In 
the right frames, the constituent contributions from df , d^ d , = Cdd,sdM> an< ^ ^sp are 
shown as the thin solid, dashed, dotted, and dash-dotted lines, respectively. The thick solid 
line is for the total EDM. 
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Figure 12: The Thallium EDM in the CPX scenario. The upper frames are for tan/5 = 5 
and the lower ones for tan/5 = 50 with $^4 = 90°. The shaded regions and lines are the 
same as in Fig. [7} 
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Figure 13: T/ie neutron EDM in the CPX scenario calculated in the CQM. The upper 
frames are for tan (3 = 5 and the lower ones for tan (3 = 50 with &a = 90° . The shaded 
regions and lines are the same as in Fig. 
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Figure 14: The neutron EDM in the CPX 
frames are for tan j3 = 5 and the lower one 
regions and lines are the same as in Fig. 
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Figure 15: The neutron EDM in the CPX scenario calculated using the QCD sum rule 
approach. The upper frames are for tan /3 — 5 and the lower ones for tan (3 = 50 with 
$a = 90°. The shaded regions and lines are the same as in Fig. [73 
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Figure 16: The Mercury EDM in the CPX scenario. The upper frames are for tan/3 = 5 
and the lower ones for tan/3 = 50 with $^4 = 90°. The shaded regions and lines are the 
same as in Fig. [721 
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Figure 17: In the upper frames, we show the Thallium EDM in the CPX scenario as 
a function of the common hierarchy factor p with several non-trivial values of ($ 1 ,$ 2 )-' 
($^$2) = (0°,0°) (solid), (90°, 0°) (dashed), (0°,90°) (dotted), and (0°,270°) (dash- 
dotted). The left frame is for tan (3 = 5 and the right one for tan f3 = 50. The lower frames 
are for the electron EDM, which makes the main contribution to the Thallium EDM, in 
the given scenario, exemplifying the case with ($1, $2) — (0°,90°) from each upper frame. 
Shown separately are the different contributions to the electron EDM from the chargino- 
(thin solid), neutralino- (thin dashed), and two-loop Higgs- (thin dash-dotted) mediated di- 
agrams. The thick solid lines are for the total EDM. We have taken $^4 = $3 = 90° in all 
frames. 



50 



Oh 

X 



10 

10 

1 



a 

u -i 
w „ 10 



Oh 



10 
10 

l 



= tan(3=5 


<*>A = *3 = 


90° 


r i ' 


<& 2 ) = (0°,90°) 




:(0°, 270°) 







I 


10 


P 




<*>A = <*>3 = 


90° 


f (0°, 270°) 

r tanp=5 


(0°,90°) 




\i («M(» ,; >. , 

(0^0°) \j 





X 




10 



w c 10 



F^A = *3 = 90° 



1 



10 




Figure 18: The neutron EDM in the CPX scenario with $^4 = $3 = 90° as a function 
of the common hierarchy factor p calculated in the CQM (upper), the PQM (middle), and 
using the QCD sum rule approach (lower) for tan /3 = 5 (left) and tan/3 = 50 (right). The 
cases with several non-trivial values of ($i,$2) ore considered as in Fig. [7?| The cases 
with (<3>i,$ 2 ) = (0°,0°) and (90°, 0°) are hardly distinguishable from each other. 
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Figure 19: Comparison of the neutron EDM calculated in the CQM (upper) and using the 
QCD sum rule approach (lower). Among the lines in Fig. [TB, the case with ($i,$2) = 
(0°,90°) is shown together with the constituent contributions: d^ d (thin solid), d^ d (thin 
dashed), d G (thin horizontal dotted), and Cbd,db (thin horizontal dash-dotted). The thick 
solid lines are for the total EDM. 
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Figure 20: The neutron EDM calculated in the PQM. In the upper frames, among the lines 
in the middle frames of Fig. [Tffi the case with (<!>!, <3> 2 ) = (0°,270°) is shown together with 
the constituent contributions: (thin solid), d® (thin dashed), and the main contribution 
from df (thin dotted). In the lower frames, the strange-quark EDM is shown as functions 
of the common hierachy factor p. The thin solid, dashed, dotted lines, and horizontal dash- 
dotted lines are for the contributions from the chargino-, neutralino- , gluino- and two-loop 
Higgs-mediated diagrams. The thick solid lines are for the total EDM. 
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Figure 21: In the upper frames, we show the Mercury EDM in the CPX scenario with 
$a = $3 = 90° as a function of the common hierarchy factor p for tan/3 = 5 (left) 
and tan/3 = 50 (right). It is hardly affected by ($i,$2), because of the dominance of the 
contribution from d^ d , see Fig. \TSi In the lower frames, we show the dominant CEDM of 
the down quark, d d , as a function of p. The lines are the same as in the lower frames of 
Fig.M 
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Figure 22: The deuteron EDM in the CPX scenario with &a = $3 = 90° as a function 
of the common hierarchy factor p for tan/3 = 5 (left) and tan/? = 50 (right). We have 
taken rf|> xp = 3 x 10 _27 ecm. It is hardly affected by ($!,$ 2 ); see the upper frames. In 
the lower frames, the case with ($!,$ 2 ) = (0°,90°) is shown together with the constituent 
contributions: d® d (thin solid), d^ d (thin dashed), C^f = Cdd,sd,bd (thin lower horizontal 
dotted), and d G (thin upper horizontal dash-dotted). The thick lines are for the total EDM. 
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Figure 23: In the upper-left frame, we show the Thallim EDM in the MCPMFV scenario 
with tan/3 = 40 as a function of $2 for several values of the CP-violating phases ($1, $3): 
(0°,0°) (solid), (0°,90°) (dashed), (90°, 0°) (dotted), and (270°, 0°) (dash- dotted). The cases 
with (0°,0°) and (0°,90°) are hardly distinguishable from each other due to the dominance 
of the electron EDM, df. In the upper-right, lower-left, and lower-right frames, we show 
the electron EDM as a function of $ 2 when ($i,$ 3 ) = (0°,90°) ; (90°, 0°), and (270°, 90°), 
respectively. The lines are the same as in Fig. [FT 
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Figure 24: In the upper- left frame, we show the neutron EDM calculated in the chial quark 
model in the MCPMFV scenario with tan j3 = 40 as a function of $2 for several values 
of the CP-violating phases ($ 1 ,$ 3 ): (0°,0°) (solid), (180°, 90°) (dashed), and (180°, 270°) 
(dotted). In the upper-right, lower-left, and lower-right frames, we show the neutron EDM 
as a function of $ 2 when ($i,$ 3 ) = (0°,0°), (180°, 90°), and (180°, 270°), respectively, 
together with its constituent contributions from d^ d (thin solid), d^ d (thin dashed), and d 
(thin dotted). The thick solid lines are for the total EDM. 
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Figure 25: In the upper-left frame, we show the neutron EDM calculated in the PQM 
in the MCPMFV scenario with tan j3 = 40 as a function of $2 for several values of the 
CP -violating phases ($i,$ 3 )- (0°,0°) (solid), (0°, 10°) (dashed), and (0°,90°) (dotted). 
All cases are not distinguishable. In this model, the neutron EDM is dominated by the 
strange-quark EDM, df. In the upper-right, lower-left, and lower-right frames, we show the 
strange-quark EDM as a function of $ 2 when ($i,$ 3 ) = (0°,0°), (0°,10°), and (0°,90°) ; 
respectively. The thin lines are for the chargino- (solid), neutralino- (dashed), gluino- 
(dotted), and Higgs-mediated (dash-dotted) diagrams, respectively. The thick solid lines are 
for the total EDM. 
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Figure 26: In the upper- left frame, we show the neutron EDM calculated using the QCD sum 
rule approach in the MCPMFV scenario with tan/3 = 40 as a function of $2 for several 
values of the CP -violating phases ($i,$ 3 ): (0°,0°) (solid), (0°,20°) (dashed), (0°,340°) 
(dotted), (0°, 270°) (dotted), and (0°, 90°) (dash- dotted). In the upper-right, lower-left, and 
lower-right frames, we show the neutron EDM as a function o/$2 when ($1,^3) = (0°,0°), 
(0°,20°) ; and (0°,90°) ; respectively. The lines are the same as in the lower frames of 
Fig. [IE 
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Figure 27: In the upper-left frame, we show the Mercury EDM calculated using the QCD 
sum rule approach in the M CP MFV scenario with tan/5 = 40 as a function of $2 for several 
values of the CP-violating phases ($i,$ 3 ): (0°,0°) (solid), (90°, 0°) (solid), (270°, 340°) 
(dashed), (90°, 10°) (dash- dotted), (270°, 350°) (dash- dotted), and (0°, 90°) (solid). In the 
upper-right, lower-left, and lower-right frames, we show the Mercury EDM as a function of 
$2 when ($i,$ 3 ) = (0°,0°), (90°, 0°), and (90°, 10°), respectively. The lines are the same 
as in the lower-right frame of Fig. QT 
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Figure 28: In the upper-left frame, we show the deuteron EDM calculated using the QCD 
sum rule techniques in the MCPMFV scenario with tan/3 = 40 as a function of $2 
for several values of the CP-violating phases (<3>i, <3>3): (0°,0°) (solid), (90°, 0°) (solid), 
(270°, 0°) (dashed), (0°, 10°) (dash- dotted), (0°,350°) (dash-dotted), and (0°,90°) (solid). 
^B XP = 3 x 10~ 27 is taken. The dominant contribution comes from the CEDM of the 
down quark, . In the upper-right, lower-left, and lower-right frames, we show the CEDM 
of the down quark as functions of $ 2 when ($i,$ 3 ) = (0°,0°), (90°, 0°), and (90°, 10°), 
respectively. The lines are the same as in Fig. \2£h 
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